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ABSTRACT: This project implements the IDENTIFICATION OF OPTIMAL DG LOCATIONS BY SINGLE DG 
PLACEMENT algorithm. This method first evaluates the voltage profile using the Newton-Raplisou me/hud and then it 
calculates the total PR loss of the system. After that by placing the DG at each bus, it evaluates the corresponding total PR 
losses and hence obtained the optimal placement ofDG for loss reduction and best suited voltage profile evaluation. 

I. INTRODUCTION 

Over the past few years, developments have been made in finding digital computer solutions for power-system load 
flows. This involves increasing the reliability and the speed of convergence of the numerical-solution techniques. The 
characteristics and performance of transmission lines can vary over wide limits mainly dependent on their system. Hence, the 
load flow method is used to maintain an acceptable voltage profile at various buses with varying power flow. The state of 
any power system can be determined using load flow analysis that calculates the power flowing through the lines of the 
system. There are different methods to determine the load flow for a particular system such as: Gauss-Seidel, Newton- 
Raphson, and the Fast-Decoupled method. 

Distributed Generation (DG) is a promising solution to many power system problems such as voltage regulation, 
power loss, etc. Distributed generation is small-scale power generation that is usually connected to or embedded in the 
distribution system. Numerous studies used different approaches to evaluate the benefits from DGs to a network in the form 
of loss reduction. 

This project implements the IDENTIFICATION OF OPTIMAL DG LOCATIONS BY SINGLE DG 
PLACEMENT algorithm. This method first evaluates the voltage profile using the Newton-Raphson method and then it 
calculates the total PR loss of the system. After that by placing the DG at each bus, it evaluates the corresponding total PR 
losses and hence obtained the optimal placement of DG for loss reduction and best suited voltage profile evaluation. 

II. POWER SYSTEM 

ll.llntroducti.on: The state of a power system and the methods of calculating this state are very important in evaluating 
the operation and control of the power system and the determination of future expansion for this system. The state of any 
power system can be determined using load flow analysis that calculates the power flowing through the lines of the system. 
There are different methods to determine the load flow for a particular system such as: Gauss-Seidel, Newton- Raphson 
Load, and the Fast-Decoupled method. 

Over the past few years, developments have been made in finding digital computer solutions for power-system load 
flows. This involves increasing the reliability and the speed of convergence of the numerical-solution techniques. In routine 
use, even few failures to give first-time convergence for physically feasible problems can be uneconomical. 

The characteristics and performance of transmission lines can vary over wide limits mainly dependent on their system. 
Hence, the load flow method is used to maintain an acceptable voltage profile at various buses with varying power flow. 

The transmission system is loop in nature having low R/X ratio. Therefore, the variables for the load-flow analysis of 
transmission systems are different from that of distribution systems which have high R/X ratio. 

II.2 Power System History: The development of the modern day electrical energy system took a few centuries. Prior to 
1800, scientists like William Gilbert, C. A. de Coulomb, Luigi Galvani, Benjamin Franklin, Alessandro Volta etc. worked on 
electric and magnetic field principles. However, none of them had any application in mind. They also probably did not 
realize that their work will lead to such an exciting engineering innovation. They were just motivated by the intellectual 
curiosity. In England, Michael Faraday worked on his induction principle between 1821 and 1831. The modern world owes a 
lot to this genius. Faraday subsequently used his induction principle to build a machine to generate voltage. Around the same 
time American engineer Joseph Henry also worked independently on the induction principle and applied his work on 
electromagnets and telegraphs. 

For about three decades between 1840 and 1870 engineers like Charles Wheatstone, Alfred Varley, Siemens 
brothers Werner and Carl etc. built primitive generators using the induction principle. It was also observed around the same 
time that when current carrying carbon electrodes were drawn apart, brilliant electric arcs were formed. The 
commercialization of arc lighting took place in the decade of 1870s. The arc lamps were used in lighthouses and streets and 
rarely indoor due to high intensity of these lights. Gas was still used for domestic lighting. It was also used for street lighting 
in many cities. 

However with the increase in load large voltage and unacceptable drops were experienced, especially at points 
that were located far away from the generating stations due to poor voltage regulation capabilities of the existing dc 
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networks. One approach was to transmit power at higher voltages while consuming it at lower voltages. This led to the 
development of the alternating current. 

As a consequence of the electric utility industry deregulation and liberalization of electricity markets as well as 
increasing demand of electric power, the amount of power exchanges between producer and consumer are increases. In this 
process, the existing transmission lines are overloaded and lead to unreliable system. The countries like India with increasing 
demand of electric power day by day it is difficult to expand the existing transmission system due to difficulties in right of 
way and cost problem in transmission network expansion. So, we need power flow controllers to increasing transmission 
capacity and controlling power flows. 

III. LOSSES IN POWER SYSTEM 
III.l Introduction: In India, average Transmission and Distribution losses, have been officially indicated as 23% of the 
electricity generated. However, as persample studies earned out by independent agencies, theselosses have been estimated to 
be as high as 50% in some states. In arecent study carried out by SBI Capital Markets for DVB, the Transmission and 
Distribution losses have been estimated as 58%. With the setting upof State Regulatory Commissions in the country, 
accurate estimation ofTransmission and DistributionLosses has gained importance as the level of losses directly affects the 
sales andpower purchase requirements and hence has a bearing on the determination ofelectricity tariff of a utility by the 
commission. 

111.2. Transmission and Distribution Losses 

III. 2.1 Losses in Transmission lines: Losses in the transmission lines can be determined less complicated compared to 
transformers and distribution systems. The basic computation of it usually surrounds to the fundamentals of ohm's law. Due 
to the simplicity of the transmission lines configuration, solving for its line losses requires no advance knowledge in any 
electrical principles. However, there are also portion of these line losses that better understanding is necessary. 
Total transmission lines losses can be broken down into three relevant parts namely; conductor losses, dielectric heating 
& radiation losses, and coupling & corona losses. It is because current flows through a transmission line and a line has a finite 
resistance there is an un-avoidable power loss. This is sometimes called conductor loss or conductor heating loss and is 
simply a power loss. 

Conductor loss depends somewhat on frequency because of a phenomenon called the skin effect. In an AC system, the 
flow of current in the cross section of the wire is not uniformly distributed. Skin effect tends to make the current flow 
concentrated more in the outer layer of the conductor. Since a very small area of the wire carries that current, line resistance 
increases at the same time increases the dissipated power. 

Corona is luminous discharge that occurs between the two conductors of a transmission line. When difference of 
potential between them exceeds the breakdown voltage of the dielectric insulator. Generally when corona occurs the 
transmission line is destroyed. If the separation between conductors in a metallic transmission line is appreciable fraction of 
wavelength. The electrostatic and electromagnetic fields that surround the conductor. Cause the line to act as if it were an 
antenna and transfer energy to any nearby conductive material. The energy radiated is called radiation loss and depends on 
dielectric material conductor spacing and length of transmission line. It reduces by properly shielding the cable e.g. STP and 
coaxial has less radiation loss It is also directly proportional to the frequency. 

7/7.2.2. Losses in Distribution lines: The term "distribution line losses" refers to the difference between the amount of 
energy delivered to the distribution system from the transmission system and the amounts of energy customers' are billed. 
Distribution line losses are comprised of two types: technical and non-technical . 

It is important to know the magnitude and causality factors for line losses because the cost of energy lost is 
recovered from customers. As a result of the composition and scale of the Hydro One distribution system it is not economic 
to provide metering and the supporting processes capable of measuring line losses directly. Since energy meters do not total 
data for the same periods, and load varies over time, a direct measurement of actual losses is not feasible. Instead, Hydro 
One relies on studies which are designed to calculate the magnitude, composition and allocation of system losses based on 
annual aggregate metering information for energy purchases, sales and system modeling methods. These studies are 
conducted with the energy assistance of industry experts in this field to ensure appropriate scientific methods and modeling 
techniques are utilized in establishing the magnitude, composition and allocation of losses. 

111.3. Components of Transmission and Distribution losses 

Energy losses occur in the process of supplying electricity to consumers due totechnical and commercial losses. The 
technical losses are due to energy dissipated in the conductors and equipment used for transmission, transformation, sub - 
transmission and distribution of power. These technical losses are inherentin a system and can be reduced to an optimum 
level. The losses can befurther sub grouped depending upon the stage of power trans formation &transmission system as 
Transmission Losses (400kV/220kV/132kV/66kV), asSub transmission losses (33kV /llkV) and Distribution losses 
(llkV/0.4kv).The commercial losses are caused by pilferage, defective meters, and errors inmeter reading and in estimating 
unmetered supply of energy. 

111.4. Reasons for Transmission and Distribution losses 

Experience in many parts of the world demonstrates that it is possible to reduce the losses in a reasonably short 
period of time and that such investments have a high internal rate of return. A clear understanding on the magnitude of 
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technicaland commercial losses is the first step in the direction of reducing T&Dlosses. This can be achieved by putting in 
place a system for accurate energyaccounting. This system is essentially a tool for energy management and helps in breaking 
down the total energy consumption into all its components. It aims at accounting for energy generated and its consumption 
by various categories ofconsumers, as well as, for energy required for meeting technical requirement ofsystem elements. It 
also helps the utility in bringing accountability and efficiency in its working. 

III.4.1 Reasons for high technical losses: The following are the major reasons for high technical losses in our country: - 

1. Inadequate investment on transmission and distribution, particularly in sub-transmission and distribution. While the 
desired investment ratiobetween generation and T&D should be 1:1, during the period 1956 -97 itdecreased to 1:0.45. 
Low investment has resulted in overloading of thedistribution system without commensurate strengthening and 
augmentation. 

2. Haphazard growths of sub-transmission and distribution system with the short-term objective of extension of power 
supply to new areas. 

3. Large scale rural electrification through long 1 lkV and LT lines. 

4. Too many stages of transformations. 

5. Improper load management. 

6. Inadequate reactive compensation 

7. Poor quality of equipment used in agricultural pumping in rural areas, cooler air conditioners and industrial loads in 
urban areas. 

7/7.4.2 Reasons for commercial losses: Theft and pilferage account for a substantial part of the high transmission and 
distribution losses in India. Theft / pilferage of energy is mainly committed bytwo categories of consumers i.e. non 
consumers and bonafide consumers. Antisocial elements avail unauthorized/unrecorded supply by hooking ortapping the 
bare conductors of L.T. feeder or tampered service wires. Some ofthe bonafide consumers willfully commit the pilferage by 
way of damaging and / or creating disturbances to measuring equipment installed at their premises. Some of the modes for 
illegal abstraction or consumption of electricity aregiven below: 

1. Making unauthorized extensions of loads, especially those having "H.P." tariff. 

2. Tampering the meter readings by mechanical jerks, placement of powerfulmagnets or disturbing the disc rotation with 
foreign matters. 

3. Stopping the meters by remote control. 

4. Willful burning of meters. 

5. Changing the sequence of terminal wiring. 

6. Bypassing the meter. 

7. Changing C. T. ratio and reducing the recording. 

8. Errors in meter reading and recording. 

9. Improper testing and calibration of meters. 

IV. LOAD FLOW STUDY 

1.1 Introduction: In a three phase ac power system active and reactive power flows from the generating station to the load 
through different networks buses and branches. The flow of active and reactive power is called power flow or load flow. 
Power flow studies provide a systematic mathematical approach for determination of various bus voltages, there phase angle 
active and reactive power flows through different branches, generators and loads under steady state condition. Power flow 
analysis is used to determine the steady state operating condition of a power system. Power flow analysis is widely used by 
power distribution professional during the planning and operation of power distribution system. 

3.2 Why we use it? 

For planning the operation of a power system, its improvement and also its future expansion require following 
studies such as load flow studies, short circuit studies and stability studies. Load flow studies are used to ensure that 
electrical power transfer from generators to consumers through the grid system is stable, reliable and economic. Load flow 
studies are most important of all power system analysis, because these are used in planning studies of power system network 
to determine if and when specific elements will become overloaded. This is important, as the magnitudes of voltages at every 
bus are required to be held within a specified limit. The objectives of any load-flow study is to provide the following 
information- 

•Voltage magnitude and phase angle at each bus. 

•Real and Reactive power flowing in each element. 
Once the bus voltage magnitudes and their angles are computed using the load flow, the real and reactive power flow 
through each line can be computed. Also based on the difference between power flow in the sending and receiving ends, the 
line losses in any particular line can also be calculated. It is helpful in determining the best location as well as optimal 
capacity of proposed generating station, substation and new lines. In the power flow problem, it is assumed that the real 
power P and reactive power Q at each Load Bus are known. For this reason, Load Buses are also known as PQ Buses. For 
Generator Buses, it is assumed that the real or active power generated P and the voltage magnitude Vis known. For the Slack 
Bus, it is assumed that the voltage magnitude V and voltage phase angle of the buses are known. In this work Newton 
Raphson method for load flow is used, because of its reliability towards convergence and not sensitive nature to the starting 
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solution. In large-scale power flow studies, the Newton raphson has proved most successful because of its strong 
convergence characteristics. 

3.3 Reactive Power: Power factor is defined as the ratio of real power to apparent power. This definition is often 
mathematically represented as Kw/Kva , where the numerator is the active (real) power and the denominator is the 
(active+reactive) or the apparent power 

Power Factor = Active power/ Apparent power = kW/kVA = Active power/ (Active Power +Reactive Power) 
= kW/(kW+kVAr) 

The higher kVAr indicates low power factor and vice versa. 

3.4 Objective of load flow study: Power flow analysis is very important in planning stages of new networks or addition to 
existing ones like adding new generator sites, meeting increase load demand and locating new transmission sites. 

• The load flow solution gives the nodal voltages and phase angles and hence the power injection at all the buses and 
power flows through interconnecting power channels. 

• It is helpful in determining the best location as well as optimal capacity of proposed generating station, substation 
and new lines. 

• It determines the voltage of the buses. The voltage level at the certain buses must be kept within the closed 
tolerances. 

• System transmission loss minimizes. 

• Economic system operation with respect to fuel cost to generate all the power needed 

• The line flows can be known. The line should not be overloaded, it means, we should not operate the close to their 
stability or thermal limits. 

3.5 Importance of maintaining bus Voltage Profile in power system: One of the major aims of system operator is to 
maintain system parameters within feasible operational margin. During certain period of time in a day-load changes very 
rapidly, especially when it tends to increase, a generalized voltage decrement and stressed conditions for reactive power 
sources occurs. Maintaining Voltage profile under varying load conditions is need of today under unpredictable load 
scenario. 



3.6 Classification of buses in load flow analysis: For load flow studies it is assumed that the loads are constant and they are 
defined by their real and reactive power consumption. It is further assumed that the generator terminal voltages are tightly 
regulated and therefore are constant. The main objective of the load flow is to find the voltage magnitude of each bus and its 
angle when the powers generated and loads are pre-specified. To facilitate this the different buses of the power system are 
classified as shown in below. 

[ 



Classification 



i ; — I ~ — ; 1 



IH.6.1 Load Buses (PQ): In these buses no generators are connected and hence the generated real power P Gi and reactive 
power Q Gi are taken as zero. The load drawn by these buses are defined by real power -P Li and reactive power -Q Lj in which 
the negative sign accommodates for the power flowing out of the bus. This is why these buses are sometimes referred to as 
P-Q bus. The objective of the load flow is to find the bus voltage magnitude IV,I and its angle <5,. 

111. 6.2 Voltage Controlled Buses (PV): These are the buses where generators are connected. Therefore the power 
generation in such buses is controlled through a prime mover while the terminal voltage is controlled through the generator 
excitation. Keeping the input power constant through turbine -governor control and keeping the bus voltage constant using 
automatic voltage regulator, we can specify constant P Gi and I V, I for these buses. This is why such buses are also referred to 
as P-V buses. It is to be noted that the reactive power supplied by the generator Q Gi depends on the system configuration and 
cannot be specified in advance. Furthermore we have to find the unknown angle <5, of the bus voltage. 

111. 6.3 Slack or Swing Bus: Usually this bus is numbered 1 for the load flow studies. This bus sets the angular reference 
for all the other buses. Since it is the angle difference between two voltage sources that dictates the real and reactive power 
flow between them, the particular angle of the slack bus is not important. However it sets the reference against which angles 
of all the other bus voltages are measured. For this reason the angle of this bus is usually chosen as 0°. Furthermore it is 
assumed that the magnitude of the voltage of this bus is known. 

Now considering a typical load flow problem in which all the load demands are known. Even if the generation 
matches the sum total of these demands exactly, the mismatch between generation and load will persist because of the 
line 1 2 R losses. Since the 1 2 R loss of a line depends on the line current which, in turn, depends on the magnitudes and angles 
of voltages of the two buses connected to the line, it is rather difficult to estimate the loss without calculating the voltages 
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and angles. For this reason a generator bus is usually chosen as the slack bus without specifying its real power. It is assumed 
that the generator connected to this bus will supply the balance of the real power required and the line losses. 

IV. Load Flow Solution 
IV.l Introduction: In Power System Engineering, the load flow study (also known as power flow study) is an important 
tool involving numerical analysis applied to a power system. Unlike traditional circuit analysis , a power flow study uses 
simplified notation such as a one line diagram and per unit system , and focuses on various forms of AC power (i.e. reactive 
,real and apparent ) rather than voltage and current. It analyses the power systems in normal steady state operation. There 
exist a number of software implementations of power flow studies. 

IV.2 Formation of y-bus matrix by Direct Inspection method: The method of building network bus admittance and bus 
impedence matrix involves transformation and inversion of matrices An alternative method for bus admittance matrix is 
based on the simple algorithm obtained by inspection of the network inter connections. 

Bus admittance is often used in power system studies. In most of the power system studies it is required to form Y- 
bus matrix of the system by considering certain power system parameters depending upon the type of analysis. 

Y-bus may be formed by inspection method, only if there is no mutual coupling between the lines. Every 
transmission line should be represented by IT -equivalent. Shunt impedances are added to diagonal element corresponding to 
the buses at which these are connected. The off diagonal elements are unaffected. The equivalent circuit of Tap changing 
transformer is included while formingY-busmatrix. 

Conseder an n bus power system (excluding the ground). taking ground as referrance and considering the single 
phase equivalent of the balanced network. 

TY n Y ln l 

Generalized Y-bus = I : : : 

Where, 

Yii=Self admittance 
Yij=Transfer admittance 

IV.2.1 Flow chart: i 1 




Compute the line admittance 

Yft = ¥iJ = R~TJx~j 



Compute the no of bus n = 
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The diagonal element Yii of the bus admittance matrix is the sum of the admittances of all the elements incident on 
bus i including the element between bus i and ground. 

The off diagonal element Yii is negative of the admittances of the element connected between buses / and / 

IV.3 Methods of load flow solution: There are mainly three methods for load flow studies: 

• Gauss Siedel method 

• Newton Raphson method 

• Fast decoupled method 

IV. 3.1 Gauss Seidel Method: The Gauss seidel method is a iterative algorithm for solving non linear algebraic equations. An 
initial solution vector is assumed,choosen from the past experience,statistical data or from practical considerations. At every 
subsequent iteration ,the solution is updated till the convergence is reached. 

IV. 3.1.1 Steps for Gauss Seidal method: Step 1: Initially assume all the buses are PQ buses,except the slack bus. This means 
that the (n-1) complex bus voltages have to be determined. 
Step 2:The Slack bus is generally numbered one. 

Step 3:PV buses are numbered in sequence and PQ buses are ordered next in sequence. 
From the equation of complex power we have, 

Si*=Vi*(EjJ =1 YikVk) 

Or, Pi-jQi=Vi*(2Xi YikVk) 

Or, ^^=(2X1 YikVk) 

Or, ^-=YiiVi+£fc=i YikVk 

Or, Vi=i * [™ - Zk=i YikVk] (a) 

Above equation is an implicit equation,since the unknown variables appears on the both sides of the equation. Hence 
it needs to be solved by an iteration technique. 

IV.3. 1.2 Algorithm For Gauss Seidel Method: Stepl: Prepare the required data. 

Step2: Formulate the Ybus matrix by the method of direct inspection. 

Step3: Assume initial voltages for all buses i= 2,3,4. ...n., (assume i=l slack bus) 

Step4: The complex bus voltagesand all n-1 buses(except the slack bus) are taken l.OzO p. u.This is normally called a flat 
start. 

Step5: Update the voltages in any (r+l)thiterations.The voltages are given by 

Vi(r+1)=^ * - Yj-Ji YikVk{r + 1) - ££=;+i YikVk(r)....{\) 

Step6: Continue the iterations till 

\AVi(r + 1)| = \Vi(r + 1) - Vl(r)\ <£ 

Where £ is the tolerance value.. Generally we use 0.0001 p.u. 
Step7: Compute slack bus power after voltages have conversed using 

Si*=Vi*(£fe =1 YikVk) 
Step8: Compute all line flows of a line (or transformer) between the buses. 



The figure represents the n model of a transmission line or a transformer. Here 
I(ik)=(Vi-Vk)y(ik)+Vi*y(iko) 
I(ki)=(Vk-Vi)y(ik)+Vk*y(iko) 
S(ik)=ViI(ik)* 
S(ki )= Vkl(ki)* 

The complex power loss in the line is given by S(ik)+S(ki). The total loss in the system is calculated by summing 
the loss over all the lines. 

IV.3. 1.3 Gauss Seidel method when p-v buses are present: Some of the buses in an n-bus power system are voltage 
controlled buses where real power(P) and bus voltage value(IVI) are specified,but reactive power Qi and voltage angle 8 are 
unknown. Hence it is necessary to first make an estimate of Qi.We can find the value of Qi by using the following 
equation- Qi= -Im{ Vi*££ =1 VikVk} 
Where ,1m stands for imaginary part. 
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At any (r+l)thiteration,at the ith PV bus 

Qi(r+1)= -lm{(VKr))*Tt=i YikVk(r + 1) + (Vi(r)) * YikVk{r)} (2)) 

The steps for ith PV bus are as follows: 

Stepl:Compute Qi(r+1) , using equation (2) 
Step2:Calculate Vi using equation (1) with Qi = Q i(r+l) 

Step3:Since IVil is specified at the PV bus, the magnitude of Vi obtained in step 2 has to be modified and set to the 
specified value IVi,spl. Therefore, 

Vi ( r+1 )= IVi ' s P'STrt 
Stated another way, 

Vi(r+l)=IVi,splz5i(r + 1) 
The voltage computation for PQ buses doesnot change. 

IV.3.1.4 Q-limit violations : If the limit of reactive power i.e. the Q limit is violated at voltage controlled bus or PV buses 
during any iteration(say,(r+l)th iteration at ith bus),the Q value is either less than the minimum Q value(Qi,min) or greater 
than the maximum Q value(Qi,max).It means that the voltage cannot be maintained at the specified value due to the lack of 
reactive power support.This bus is then treated as a PQ bus in the (r+l)th iteration and the voltage is calculated with the 
value of Qi set as follows: 

If Qi<Qi,min then Qi=Qi,min 

If Qi>Qi,max then Qi=Qi,max 
If in the subsequent iteration, Qi falls within the limits, the bus can be switched back to PV status. 

IV.3.1.5 Acceleration of convergence: The number of iterations increases with the increase in the size of the system, while 
we use Gauss Seidel method. The number of iterations required can be reduced if the correction in voltage at each bus is 
accelerated, by multiplying with a constant a.cris called the acceleration factor. 

In the (r+l)th iteration we can write, 

Accelerated Vi(r+1) = Vi(r) +a(Vi(r+l)-Vi(r) ), 

Where, a is a real number. When a=l, the value of Vi(r+1) is the computed value. If 1< a <2, then the value 
computed is extrapolated. Generally, a is taken between 1.6 to 2. OA wrong value may lead to divergence of the solution. 

At PQ buses (pure load buses) if the voltage magnitude violates the limit. It simply means that the specified reactive 
power demand cannot be supllied, with the voltage maintained within acceptable limits. 
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Let Qi(r+1 


) =Qi, max 




LetQi(r+l)=Qi, min 
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Compute Vi (r+1) 




Compute 

5i(r+l)=ZVi(r+l) 

Vi(r+l)=IV,isplz5i(r+l) 



Calculate line flows, 
slack bus voltage, losses 
and display 



IV.3.2 Newton Raphson Method: In application of the NR method, we have to first bring the equations to be solved to the 
form f(xl,x2,...xn)=0, where xl,x2..xn are the unknown variables to be determined. Let us assume that the power system has 
nl PV buses and n2 PQ buses. In polar co ordinates the unknown variable to be determined are. 



bus i,at all the PV and PQ 



IV.3.2. 1 Steps for Newton-Raphson method: Stepl: Si the angle of the complex bus voltage a 
buses. This gives us nl+n2 unknown variables to be determined. 

Step2: IVil, the voltage magnitude of bus i, at all PQ buses. This gives us n2 unknown variables to be determined. 

Therefore the total number of unknown variables to be computed is nl+2n2 for which we nl+2n2 consistent 
equations to be solved. The equations to be solved are given by 

APi = Pi,sp - Pi, cat = 0 

AQi = Qi,sp - Qi.cal = 0 

Where, 

AP = Active power residue 

AQ= reactive power residue 
Pi, sp =specified active power at bus i 
Qi, sp =Specified reactive power at bus i 
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Pi, cai=calculated value of active power using voltage estimates 
Qi, caZ=calculated value of reactive power using voltage estimates 

IV. 3.2.2 Algorithm for N-R method: Stepl: Formulate Bus Admittance Matrix Y bus. 
Step2: Assume initial voltages as follows 

Vi=IVi, spl^O (at all PV Buses) where Vi,sp=Specified voltage value 
Vi=lzO (at all PQ buses) 

Step3: At (r+l)th iteration calculate Pi(r+l),at all the PV and PQ Buses and calculate Qi(r+1) at all the PQ buses using Vi(r) 
values. 

The formulae to be used are 

Pi,cal=GiilVt| 2 + Yl=i\Vi\\Vk\{Gik * cos Sik + Bik * sinSik) 

Qi,cal= BiilVi | 2 + =i I Vi \ \ Vk \ (Gik * sin Sik -Bik* cosSik) 

Step4: Calculate the power mismatches i.e. Power residues 
APi(r) = Pi,sp - Pi,cal(r + 1) (at PV and PQ buses) 
AQi(r) = Qi,sp - Qi,cal(r + 1) (at PQ buses) 
Step5: Calculate the Jacobian [J(r)] using Vi(r) 

A5(r) Ap(r) 
Step6: Compute ajv^i = ^* A(?(r) 

Step7:Update the variables as follows- 

Si (r + 1) = Si (r) + ASi (r) (at all buses) 
IVi(r+l)l=IVil(r)+A|Fi|(r) 
Step8: Go to step 3 and iterate till the the power mismatches are within the acceptable limits. 



3.3.2.3 Formation of Jacobian: 



















A.Q" 










<5X3„ 




*K 1° 



There are some formulas given below for making the Jacobian in NR load flow method, theseare... 

D^=-Qi-Bii|^| 2 

2) g= \Vi\\Vk\(Gik * sin Sik - Bik * cosSik) 

3) ^\Vk\= \Vi\\Vk\(Gik * cosSik + Bik * sinSik) 

4) ^y|7i| =Pi+Gii|Kt| 2 

5) f=Pi-Gii|^| 2 

6) |J = -\Vi\\Vk\(Gik * cosSik + Bik * sinSik) 

7) ^|7i|=Qi-Bii|7i| 2 

8) ^y|Kfc| = \Vi\\Vk\(Gik * sinSik - Bik * cosSik) 
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IV.3.2.4 Flow Chart for Newton Raphson method 



Form bus admittance 



Assume initial bus 

voltage Vi ; 
i= 1, 2 ...n, \ + slack 



Compute Pi, caland Qi, cal 
i = 2, ,n, r/slack 



Set iteration count r =0 
tolerance, €=0.0001 



Compute power residues 
APi,r = Pi,sp-Pi,calAQi,r 
= Qi ,sp-Qi, cal 



Calculate elements 
of Jacobian 



Calculate change in 
angle and voltage 



Update variables 5(r+l) 
=5(r) + A 5(r) 
|V|(r+l)=|V|(r)+AV(r) 




Calculate line flows, 
slack bus voltage, 
loses and display 
result 



No 



Go to A 



IV. 3.3 Fast Decoupled method 

If the co efficient matrices are constant,the need to update the jacobian at every iteration is eliminated. This has 
resulted in development of fast decoupled load flow. Assumptions used in this method are- 

• Bij » G t j (Since,the X/R ratio of transmission lines is high in well desighned system. 

• Real power changes are less sensitive to voltage magnitude changes and are most sensitive to changes in phase 

angle A5. 
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i.ecos(5[ — Sj)~l and sin(5j — 8j)~0 

• Similarly, reactive power changes are less sensitive to changes in angle and are mainly dependent on changes in voltage 
magnitude. 

• Therefore the element of theJacobian matrix become 

H ik= L ik = -\Vi\\Vk\Bik(i * k) 
H u = L u = -Bii\Vi\ 2 

With this value we can perfom the load flow analysis by taking less time. 

V. DISTRIBUTED GENERATION 
V.l Introduction to Distributed Generation 

What is Distributed Generation? 

Distributed generation is an approach that employs small-scale technologies to produce electricity close to the end 
users of power. DG technologies often consist of modular (and sometimes renewable -energy) generators, and they offer a 
number of potential benefits. In many cases, distributed generators can provide lower-cost electricity and higher power 
reliability and security with fewer environmental consequences than can traditional power generators. Distributed generation, 
also called on-site generation, dispersed generation, embedded generation, decentralized generation, decentralized energy or 
distributed energy, generates electricity from many small energy sources. Distributed generation reduces the amount of 
energy lost in transmitting electricity because the electricity is generated very near where it is used, perhaps even in the same 
building. This also reduces the size and number of power lines that must be constructed. 

In contrast to the use of a few large-scale generating stations located far from load centers— the approach used in the 
traditional electric power paradigm— DG systems employ numerous, but small plants and can provide power onsite with little 
reliance on the distribution and transmission grid. DG technologies yield power in capacities that range from a fraction of a 
kilowatt [kW] to about 100 megawatts [MW], Utility-scale generation units generate power in capacities that often reach 
beyond 1,000 MW. 

Distributed generation takes place on two-levels: the local level and the end-point level. Local level power 
generation plants often include renewable energy technologies that are site specific, such as wind turbines, geothermal 
energy production, solar systems (photovoltaic and combustion), and some hydro -thermal plants. These plants tend to be 
smaller and less centralized than the traditional model plants. They also are frequently more energy and cost efficient and 
more reliable. Since these local level DG producers often take into account the local context, the usually produce less 
environmentally damaging or disrupting energy than the larger central model plants. 

V.2 Benefits of Distributed Generation 

What are the Potential Benefits of DG Systems? 

Consumer advocates who favor DG point out that distributed resources can improve the efficiency of providing 
electric power. They often highlight that transmission of electricity from a power plant to a typical user wastes roughly 4.2 
to 8.9 percent of the electricity as a consequence of aging transmission equipment, inconsistent enforcement of reliability 
guidelines, and growing congestion. At the same time, customers often suffer from poor power quality — variations in voltage 
or electrical flow — that results from a variety of factors, including poor switching operations in the network, voltage dips, 
interruptions, transients, and network disturbances from loads. Overall, DG proponents highlight the inefficiency of the 
existing large-scale electrical transmission and distribution network. Moreover, because customers' electricity bills include 
the cost of this vast transmission grid, the use of on-site power equipment can conceivably provide consumers with 
affordable power at a higher level of quality. In addition, residents and businesses that generate power locally have the 
potential to sell surplus power to the grid, which can yield significant income during times of peak demand. 

Beyond efficiency, DG technologies may provide benefits in the form of more reliable power for industries that 
require uninterrupted service. The Electric Power Research Institute reported that power outages and quality disturbances 
cost American businesses $119 billion per year. In 2001, the International Energy Agency (2002) estimated that the average 
cost of a one-hour power outage was $6,480,000 for brokerage operations and $2,580,000 for credit card operations. The 
figures grow more impressively for the semiconductor industry, where a two hour power outage can cost close to 
$48,000,000. Given these numbers, it remains no mystery why several firms have already installed DG facilities to ensure 
consistent power supplies. 

Perhaps incongruously, DG facilities offer potential advantages for improving the transmission of power. Because 
they produce power locally for users, they aid the entire grid by reducing demand during peak times and by minimizing 
congestion of power on the network, one of the causes of the 2003 blackout. And by building large numbers of localized 
power generation facilities rather than a few large-scale power plants located distantly from load centers, DG can contribute 
to deferring transmission upgrades and expansions — at a time when investment in such facilities remains constrained. 
Perhaps most important in the post-September 1 1 era, DG technologies may improve the security of the grid. 

Environmentalists and academics suggest that DG technologies can provide ancillary benefits to society. Large, 
centralized power plants emit significant amounts of carbon monoxide, sulfur oxides, particulate matter, hydrocarbons, and 
nitrogen oxides. 

Finally, DG can help the nation increase its diversity of energy sources. Some of the DG technologies, such as 
wind turbines, solar photovoltaic panels, and hydroelectric turbines, consume no fossil fuels, while others, such as fuel cells, 
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microturbines, and some internal combustion units burn natural gas, much of which is produced in the United States. The 
increasing diversity helps insulate the economy from price shocks, interruptions, and fuel shortages 

V.3 Usage of Distributed Generation 

How Much are Renewable and DG Systems Used in the American Electric Utility Sector? 

Despite the immense environmental, technical, and financial promise of renewable energy systems, such generators 
still constitute a very small percentage of electricity generation capacity in the United States. Throughout the 1970s, some 
policy experts expected renewable energy systems to be used for much more generation capacity than they have. Dr. Arthur 
Rosenfeld, one of the five CEC commissioners serving from 2002 until the present, noted that President Carter had told him 
(during his presidency in the late 1970s) that he expected renewable energy systems to reach 10 percent of national 
electricity capacity by 1985. However, Carter's expectation went unfulfilled: excluding large hydroelectric generators, 
renewable energy technologies in 2003 comprised only about 2 percent of the U.S. electricity generation mix 

The relatively minor use of renewable energy systems has created a general attitude among energy analysts, 
scholars, and laboratory directors that the technologies are not viable sources of electricity supply. For example, 
RodeySobin, former Innovative Technology Manager for the Virginia Department of Environmental Quality, argues that "in 
many ways, renewable energy systems were the technology of the future, and today they still are." Ralph D. Badinelli, a 
professor of Business Information Technology at Virginia Tech, explains that renewable energy technologies do not 
contribute significantly to U.S. generation capacity because "such sources have not yet proven themselves ... Until they do, 
they will be considered scientific experiments as opposed to new technologies." Similarly, Mark Levine, the Environmental 
Energy Technologies Division Director at the Lawrence Berkeley National Laboratory, comments that despite all of the hype 
surrounding renewable energy, such systems are still only "excellent for niche applications, but the niches aren't large." 

DG/CHP technologies have an only slightly better record. In 2004, the Energy Information Administration 
characterized only 3.1 percent of electricity generation capacity as commercial or industrial combined heat and power 
(33,217 MW out of 1.49 terrawatts [TW]). The EIA also estimated that in 2002 only 0.9 gigawatts (GW) of distributed 
generation capacity existed in the United States. Similarly, the EIA's 2005 Annual Energy Outlook projected that CHP 
systems are not widely used in the electric power sector, amounting to 0.053% of utility generation (197 billion kWh out of 
3,700 billion kWh). Tom Casten, the Chair and Chief Executive Officer of Primary Energy, a manufacturer of fuel 
processing cogeneration steam plants, notes that even though CHP plants can reduce energy costs for industrial firms by over 
40 percent, such plants remain "the exception instead of the rule." 

V.4 Distriduted Generation in our system 

5.4.1 Algorithm: The computational steps involved in finding the optimal Distribution Generator (DG) size and location 
to minimize the loss in a radial distribution system are summarized in following: 

1. Perform the load flow analysis by using Newton Raphson method for determining the voltage profile and Total Loss for 
Radial Distribution system. 

2. In load flow analysis obtain the branch current 7 i; between two buses by using , 

/y = -j- £■ Where, V t — Voltage of bus i 

Vj = Voltage of bus j 
R tj = Resistance between bus i and j 

3. Obtain the active component, I a and reactive component, I r of the branch currents /y . 

4. We calculate total I%R loss by using, 

Where, P La = Yu=i ^ai ^ij Due to active component of the current, 
P Lr = Yi=i & Rij Due to reactive component of the current, 

5. We calculate the power loss P™ w associated with the active component of branch current when DG is connected.lt is 
given by 

P La W =Yi=iUai + DiI DG ) 2 R t where ZJ £ =1; if branch iea 
=0,otherwise 

' a' is the set of branches connected between the source and the bus m where DG is placed. 

6. Repeat steps 1 to steps 5 and calculate P ia (new) by placing DG at each bus. 

7. Calculate the power saving by applying the formula given below 

S=P La -P™ w 

8. The DG current IDG that provides maximum saving can be obtained from 
-22r=i(A4i + D i I DG )R i =0 

9. Calculate DG current for maximum power saving and is given by 

, _ Z iea laiRi 

10. Calculate the corresponding DG Size to be placed at each bus ,is given by 
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At Bus m, 

Pdg=KiIdg > where V m is the Voltage Magnitude of Bus 'm' 

11. Update the active power component 

Pnew = PgEN + ^DG ~ PiOAD 

12. Perform load flow analysis with updated active power component obtained in step 11. 

13. Perform the step 10 to step 12 , and record the Voltage Profiles and Total line loss. 

14. Obtain the Optomal location of Distributed Generation (DG) for Total line loss minimization and Voltage profile 
improvement. 

V.4.2 Flow Chart 

| Start ~j 



Read System Data 
P, Q, Vsp, 5, R, X, count 

Jl 



Execute Newton Raphson 




Record Voltage profile & 
Total I 2 R line loss Without DG 



II ~ 

Set count=2 
(place the DG at bus no. 2) 

it 

Calculate Branch Current I DG from the source to the 
particular BUS 

Calculate DG power 

P D G ^ count X he 



Update new active Power 

Pnew = PGEN + P DG ~ P L0AD 




Execute load flow (NR method) with updated 
P and Obtain new Voltage Profile and I 2 R 
Line Losses. 



I 

Compute the size of DG for 
Multiple locations 

II 

Record Voltage profile & 
Total I R line loss 
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I . 

Display Voltage Profiles & Total I R line 
loss with & without DG 

n 

Stop 



VI. TEST SYSTEM AND RESULTS 
VI. 1 Test System (12-BUS Radial System) 





Figu 



■ .4. Jinn's si JD£sz--£c:i: 



Number of buses in the system 

Number of slack bus = 1 

Number of PV buses =11 

Number of PQ buses = 0 

Number of lines in the system =11 

Maximum convergence value = 0.0001 
Base value = 1000 
Coding used for Buses: 

Slack Bus - 1 ; P V Bus - 2; PQ Bus - 3 

VI.1.1 SPECIFICATION OF TEST SYSTEM: 



Bus 


Type 


Vsp (in pu) 


Power angle 


PGi (in kW) 


QGi 

(in kVAR) 


PLi 

(in kW) 


QLi 
(in kVAR) 


1 


1 


1.06 


0 


0 


0 


0 


0 


2 


3 




0 


0 


0 


60 


60 


3 


3 




0 


0 


0 


40 


30 


4 


3 




0 


0 


0 


55 


55 


5 


3 




0 


0 


0 


30 


30 


6 


3 




0 


0 


0 


20 


15 


7 


3 




0 


0 


0 


55 


55 


8 


3 




0 


0 


0 


45 


45 


9 


3 




0 


0 


0 


40 


40 


10 


3 




0 


0 


0 


35 


30 


11 


3 




0 


0 


0 


40 


30 


12 


3 




0 


0 


0 


15 


15 



VI.1.2 LINE DATA: 



From Bus 


To Bus 


Resistence, R (in Q) 


Reactance, X (infl) 


1 


2 


1.093 


0.455 


2 


3 


1.184 


0.494 


3 


4 


2.095 


0.873 
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VII. RESULTS & DISCUSSIONS 

Performing the experiment in MATLAB Environment for Optimal Placement of DG with loss minimization technique 
and Voltage Profile improvement, we are having the results as described below: 



1. Optimal Placement of DG for Loss Minimization: 

For finding the optimal location of DG for Total Loss Minimization we placed single DG at each Bus and Total loss 
obtained in each step is compared as such. 





Total Loss with DG placement 






















1 
















No 2 3 4 5 6 7 8 9 10 11 12 
DG 

BUS 



Fig : Optimal DG location for Total Line Loss minimization 

From the graph we observe that when we place DG at BUS No. 4 the Total Line Loss compared to DG placement 
at other buses is minimum. The Elapsed time for execution of the MATLAB program is 0.385176 seconds with 
total line loss = 2.91113^10- 6 p.u. 

So, we can suggest that for Total Line Loss reduction in our test system, the optimal location for the single DG 
placement is at BUS No. 4 



2. Optimal Placement of Voltage Profile Improvement: 

For finding the optimal location of DG for Voltage Profile improvement we placed single DG at each Bus and 
compared the obtained voltage profile with the reference value (1 p. u) 



i S 9 i bus io r 



00413 1.00415 



,00505 1.00501 



- "-.:.r 



,01272 1.01263 



,02102 1.O2097 



,04275 1.Q4272 



,04343 | 1.04341 



Table : Voltage Profile before and after placing DG 
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COMPARISION OF VOLTAGE PROFILES 




[1]. In the graph Blue bar represents the voltage profile of the 12 Bus radial distribution system without DG placement and Red bar 

represents the voltage profile after placing the DG at bus no. 3 
[2]. Comparing the two Voltage Profiles with reference value (1 pu), we observed that , when the DG is placed at Bus No. 3 the 

obtained voltage is much nearer to the reference value. 
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